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, We analyze the radiative capture of the negative muon in hydrogen using 

amphtudes derived within the chiral Lagrangian approach. Besides the lead- 
ing and next to leading order terms, given by the well-known Rood-Tolhoek 
Hamiltonian, we extract from these amplitudes the corrections of the next 
order in 1/M (M is the nucleon mass). In addition, we estimate within the 
^ . same formalism also the A (1232) isobar excitation effects and processes de- 

I scribed by an anomalous Lagrangian. Using of the parameters of the model 

obtained from the analysis of the pion photoproduction, which restricts the 
arbitrariness of the vrA^A and 7A^A vertices off-shell, allows us to explain 
two times more of the discrepancy between the value gp*^^*^' of the induced 
pseudoscalar form factor gp, predicted by partial conservation of the axial 
current, and of gp extracted from the recent TRIUMF experiment, than the 
standard approach. Varying these parameters independently, one can remove 
the discrepancy completely. 
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I. INTRODUCTION 



It is well known [|T| that the charged weak interaction of the nucleon with a lepton is described 
by the weak hadron current 



where the vector part is given by the matrix element of the isovector Lorentz 4-vector current 
operator between the nucleon states, 



Jv,M) = 'l9vi<lih,-^-^^,-^^i^] (1-2) 



and the axial-vector part is analogously, 

l^..^t; -I- 1— 11 u /o 

mi 12 



JIM) = ' (-9Aiqlh,l, + ^^4^^1.75 1 ^ . (1.3) 



Here a is the isospin index, mi is the lepton mass and the 4-momentum transfer is given by 
Qiii = p'f_i — Pii, where p'^ (p^) is the 4-momentum of the final (initial) nucleon. 

The least known of the four form factors entering the currents Eqs. (|1.2| ) and (|1.3| ) is the 
induced pseudoscalar form factor, gp{qi), in the axial-vector current Jl^^- Actually, its pres- 
ence in the axial-vector current ( |1.3D tests our understanding of the basic strong and weak 
interaction processes, such as the strong ttNN vertex and the weak pion decay. Elementary 
calculations lead to 

gp{q^) = -2g^NNUmi Al{q^) , (1.4) 



where A^(gf) is the pion propagator, g^NN = 13.05 is the pseudoscalar vrA^A^ coupling 
constant and = 92.4 MeV is the pion decay constant. 

The matrix element of the axial current ^ should satisfy partial conservation of the axial 
current (PCAC). It is easy to obtain that 

75 Y^i(p) • (1-5) 



u{p'] 



u{p') 



2MgAFA{ql) - 



gpiq? 



mi 



It is seen from this equation, that if 



9p{(1i) = -2gnNNf7Tmi— 1 + 



MgA 

gnNNfm 



FA{ql 



:i.6) 
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is subtracted from gp{qi), then indeed, PC AC is valid. Here we put 

9A{qi) = 9AFA{qi) , 



(1.7) 



with = 5'a(0) = —1.267. In the chiral model @J^, the axial form factor is of the monopole 
form 



Fa{Qi 



m 



ax 



9 I 2 ' 



:i.8) 



where m^^ is the mass of the axial vector meson ai(1260). Then for the gp{qi ), Eq. ([TT 
we have 



9p{(1i) = -'^gnNNUmi-^ 

Qi 



9i + 1 + r 

\ 9ttNNJtt, 



m„ 



:i.9) 



However, this equation for gp{qi) cannot be used, because of the singularity for = 0, which 
shows its presence in the hadron radiative part of the radiative muon capture amplitude 
(RMC) for large photon momentums k and it is close to the physical region for the ordinary 
muon capture (OMC), because of the large value of the axial meson mass. So our model 
cannot be used beyond the exact Goldberger-Treiman relation and we take 



9p{qi 



(1.10) 



which is in agreement with 

The best way to search for the effect of the form factor gp{qi) is the muon capture. In the 
elementary process of OMC in the hydrogen. 



H + p — > i^f, + n, 



;i.ii) 



according to Eq. ( |1.4|) , the value of the induced pseudoscalar form factor gp is 



gT^'^V) ^ gp{ql = 0.877mJ) 



0.877m2 + ml 



6.87 c/A = -8.71 



;i.l2) 



and for 5'p(gf), Eq. ( |1.10|) , we have 



g$''^ip) = ~gp{Q[' = 0.877m^ 



"^gnNNfiT 



0.877m^ + m^^ 



0.13 c/A = -0.16^ 



(1.13) 



which is a correction of ~ 2% to the gp^'"{p), Eq. ( |1.12| ). The resulting value is 



(1.14) 



The axial form factor of the nucleon has recently been measured by the p{e, e'7c^)n reaction 
in Ref. P). The dipole form of the form factor was used and the extracted axial mass 
rriA = 1.077 ± 0.039 GeV. For this form factor, 



'2i7Th^ _|_ Q 2 

gpiQi) = -'^gnNNfTrm-r^, — t^- (1-15) 



and analogously with Eqs. ( |1.13| ) and ( |1.14| ) we have 

-g^^'^ip) = OMgA = -0.43, (1.16) 

and 

g^^^^ip) = -8.28, (1.17) 

respectively. Both values of gp'"^'"{p) are in reasonable agreement with the calculations of 
gPCAc^p\^ within the framework of the heavy baryon chiral perturbation theory (HBChPT), 



Very flat dependence of the capture rate on gp in the OMC by the proton, Eq. ( |LT1| ) 



provides its world average value with an error of ~ 20 % and particular experiments have 
an error larger by a factor of ~ 2. 



Recently, a very precise experimental study of the muon capture by ^He [1C,11], 

^^- + ^He — > v^ + '^H, (1.18) 

yields the transition rate 

Te^p = 1494 ± 4s^S (1.19) 



which allowed |T^ an extraction of the value oi gp with an accuracy of ~ 20 % from this 



experiment alone. 



1.05 ±0.19, (1.20) 



where for the reaction Eq. ( [1.18| ) 

g^^^CHe) ^ g°/'^{ql = 0.954 m^) = 6.68^^. (1.21) 



Analogously with Eq. ( |1.14|) , this value oi gp differs slightly from the one demanded by 
PCAC by about the same amount as for the reaction Eq. ( 1.11 ] 



A contribution of 20 



due to the meson exchange current effect turned out to be essential to get the calculated 
transition rate 



th 



1502 ± 32 s" 



:i.22) 



into agreement with the data Eq. (|1.19|) . Let us note that a further improvement of the 
extracted value of gp is hindered by an uncertainty of ~ 2 % in calculations [|T^ which will 
be difficult to improve. The main uncertainty arises from the less known parameters of the 
A excitation processes. 



Another attractive tool to extract the value of gp'"^'" is the RMC by the proton, 

H~ + p — > z/^ + 7 + n . 



;i.23) 



As is well known [|^, the RMC amplitude contains the pseudoscalar form factor gp in the 



form 



—X 



xgr"" 



0.877ml + ml 



1 - 



(1.24) 



where we implemented the correction of Eq. (|1.6| ) and 



L{N)\2 



l/((g-(^))^ + m^^ 



:i.25) 



for the chiral model |0J^ and 



(1.26) 



for the dipole form of the form factor. As it is seen from Eq. ( p..24| ), the form factor gp 



depends either on the square of the four-momentum transfer q 



p—p 



u + k — jj 



V — jjL = q —k (for the hadron 



-rN\2 



m^,, which 



2 

3. This enhancement makes the 



(characterizing the muon radiation process) or on the 
radiation). For large photon momentums k, {q^Y ~ + 
enhances the hadron radiation amplitude by a factor of 
reaction (|1.23| ) particularly interesting. On the other hand, the dependence on gp of the 
effective form factors gi, entering the effective RMC Hamiltonian, appears only from the 
order (9(1/M), which makes the isolation of the dependence of the photon spectrum on gp 
difficult^. The factor x in Eq. ( p..24| ) is used to study the change of the photon spectrum and 
capture rates by scaling gp. 



^In Ref. it is proposed to isolate the effect due to the hadron radiative amplitude in a very 
difficult polarization experiment. 
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The theory of the RMC was elaborated by many authors during the decades (see Refs. 0, 
T3|,15-^ and references therein). 



The nuclear Hamiltonian, suitable for use in nuclear physics calculations of the RMC pro- 



cesses, was provided first by Rood and Tolhoek [|T^]. It contains the leading and the next to 
leading order terms in 1/M derived from the conserved RMC amplitude given by a set of 
the Feynman diagrams. Christillin and Servadio rederived in an elegant way the RMC 
amplitude obtained earlier by Adler and Dothan using the low energy theorems. This 
amplitude is written in terms of elastic weak form factors and pion photoproduction ampli- 
tude, up to terms linear in k and q. It was also found |^ that higher order terms cannot 
be obtained using this method. Recently, this amplitude was produced ^ from a chiral 
Lagrangian of the Nnpuiai system. It satisfies the corresponding continuity equations and 
the consistency condition exactly. Higher order terms follow without any restriction. It was 
shown that the leading order terms coincide with those given by the low energy theorems. 
However, higher order terms differ, which is given by a different prescription to pass towards 
higher energies. 



The above mentioned set of the relativistic Feynman diagrams was used by Fearing to 



calculate the photon energy spectrum for the reaction Eq. (|1.23|) . This work was later ex- 
tended by Beder and Fearing p6|] by considering also the contribution from the A excitation 
processes. A recent comparison of the TRIUMF experiment [^,0 with the Beder-Fearing 
calculations provided a value of g'p^'"{p) which is enhanced by ~ 50% in comparison with 
the value of Eq. ( p..l2| ), which corresponds to using gp from Eq. ( |1.24| ) with x = 1.5 . This 
is the so-called ^gp puzzle'. 



In connection with the presence of the factor x in Eq. ( |1.24|) , it should be noted that 



(i) Referring only to the change of gp^^'"{p) is confusing. As it is seen from Eq. (|1.24|) , the 
whole form factor gp is scaled. 

(ii) In the experiment |2^j2^, the high energy part of the photon spectrum is measured. 
Then increasing x simulates processes enhancing this part of the spectrum[[ 

In our opinion, the variation of x can be considered as a tool to study an uncertainty in our 
knowledge of gp due to a restricted experimental accuracy. Any real difference from x = 1 
would mean violation of PCAC. 

Searching for the processes enhancing the high energy part of the photon spectrum has re- 
cently been performed within the concept of HBChPT by several authors [p9|-p3[]. Ando and 



Min considered one-loop order correlations to the tree approximation and confirmed the 



existing discrepancy. Bernard, Hemmert and MeiBner (BHM) calculated [03 both ordinary 



^ See also the discussion in Ref. |29|. 



6 



and radiative muon capture on the proton in an effective field theory of pions, nucleons and 
A isobars by using the small scale expansion |34|. According to |2^, the most probable 
explanation of the problem is a combination of many small effects. Besides the photon 
spectra, BHM present the numerical results also for the singlet (Ag) and triplet (A^) capture 
rates. This will enable us to compare our calculations with those by BHM in a more detail. 
Here we only note the difference of ~ 10% in A^. As we shall see later, about half of this 
difference arises from the use of an approximate equation for the neutrino energy in Ref. 



In Ref. |22[ , a possible explanation of the discrepancy was suggested that a fraction of spin 
3/2 orthomolecular p/ip state in liquid hydrogen can exist. The analysis of the experimental 
photon spectrum p7| , |28[] yielded 10 to 20% of this state. However, this is in sharp contrast 
with the existing calculations [^,^, which give zero fraction of this state. As noted very 
recently in Ref. [§^0, a new analysis restricts the fraction of spin 3/2 orthomolecular p/xp 
state to at most 



70. 



Finally, let us comment on Ref. ||3^ , where Cheon and Cheoun reported on the derivation of 
an additional term from a chiral model, which does not appear in the standard approach to 
the RMC on the proton and which generates a large contribution to the photon spectrum. 
As it was shown in Ref. the report ||3^ suffers from two flaws. First, the derivation 
of this terms contains an algebraic error due to an incorrect application of the covariant 
derivative in Eq. (14). After removing it, the effect is reduced by a factor of ~ 5. The 
second flaw in Ref. [|3^ is related to the introduction of the pseudovector nNN coupling by 
the vertex Ci of Eq. (18), which yields the desired term. However, the equivalent passage 
from one type of the ttNN coupling to another one is guaranteed only by the Foldy-Dyson 
unitary transformation. As shown in Ref. ||38|, when this transformation is applied to a chiral 
model with the pseudoscalar ttNN coupling, the pseudovector tcNN coupling appears in the 
resulting Lagrangian, which does not contain the incriminating term, however. Besides, the 
presence of this term in the RMC amplitude violates the Ward-Takahashi identity derived 
in Ref. . Later attempt to improve the situation |3^ suffers from the same shortcomings. 
The report 



was also criticised in Ref. |^ 



This situation makes the expectation of the result from the next TRIUMF experiment in 
helium, 



;i.27) 



with a particular tension. However, one should keep in mind also complications analogous 
to those in reaction Eq. ( |1.18| ) and non-negligible meson exchange current effects are to 
be expected, which makes the analysis much more difficult. It is clear that the Beder- 
Fearing relativistic formalism is not applicable in calculations with the realistic 3N wave 



^See also a discussion in Ref. 
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functions and a consistent non-relativistic approach should be developed. Here we make an 
independent step in this direction by making the non-relativistic reduction of the amplitudes 
derived in Ref. P] from a chiral Lagrangian of the Nnpuai system. As a result, we get an 



effective Hamiltonian which is close to the one obtained by Rood and Tolhoek |]19| but not 
identical with it. We also apply the constructed effective Hamiltonian to compute both 
the capture rates and the photon energy spectra for the reaction Eq. ( |1.23|) and for various 
spin states. Our reduction provides more terms of the order (1/M) and (1/M^) than Rood 
and Tolhoek present. Added to the leading order terms, they should reproduce with a 
good accuracy results given in Ref. Another set of terms of the order (9(l/m^) ^ 

(9(1/M^) is produced by reduction of additional relativistic amplitudes following from our 
chiral Lagrangian. We shall call it hard pion (hp) correction. Numerically, it enhances the 
photon spectra by 2-4%. 



Next we include the A isobar using again the formalism of chiral Lagrangians developed 
in Refs. 0,0, which we extend by adopting results of Refs. p3|-|45[]. Then the resulting 
NAT^pai Lagrangian consists of three terms and is characterised by three couplings and four 
arbitrary parameters A, X, Y, Z. In its turn, each term contains a tensor of the form 



;i + 4B)A + B 



1fj.1v 



B = X,Y, Z, 



[1.28) 



which ensures the independence of the A contribution to the S-matrix on the parameter 
A. The choice A = —1 simplifies the A propagator considerably. The parameters X, Y, Z, 
which reflect the off-shell ambiguity of the massive spin 3/2 field were found |]^3|-^ by 
analyzing the data on pion photoproduction]^. The value of these parameters depends on 
how the pion photoproduction amplitude is unitarized. This model does not require the use 
of the Breit-Wigner form of the A propagator. 



In the calculations of the A excitation effect in the reaction Eq. ( |L23| ), Beder and Fearing 
26| took a model for needed vertices with Q^^, = 6fj_i,, the Breit-Wigner form of the A 
propagator and the needed 7A^A coupling from Ref. MM , thus introducing an inconsistency 



into calculations. This model provides about 7% effect from the A excitation to the photon 
spectrum. 

The vriVA and 7A^A vertices including the off-shell parameters X, Y and Z were discussed 
in Ref. and the vrA^A vertex of the form Eq. (|1.28|) was considered also in the small 
scale expansion [Q. However, the dependence of the results on these parameters was not 
exploited in any of the calculations performed within the framework of HBChPT. 



Let us note that besides the adopted model [^^-^, other models |pD|P| were developed 



^These corrections up to the order (1/M^) were discussed in Ref. [41|. 

^This model describes well also the latest data on the vr*^ electroproduction on the proton [48| 
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to describe the production of pions on protons by the electromagnetic interaction. All 
these models consider the same non-resonant Lagrangian of the Nirpu system, but differ 
principally in the treatment of the A isobar and in the method of unitarization of the vrA^ 
amplitude. 

We have also analyzed the contribution due to amplitudes constructed from an anomalous 
Lagrangian of the npujai system We have found that the influence of this contri- 

bution on the photon energy spectrum is not significant. Earlier estimate of a contribution 
which arises from the Wess-Zumino-Witten part of the anomalous Lagrangian was reported 
in Ref. 



One can find in the literature an attempt to study the form factor gp the reaction of 



electroproduction of charged soft pions off the proton [p^,p5 



e + p 



e + vr + n , 



;L29) 



The starting point of this attempt is the soft pion production amplitude given as 
hMl\q,k) '-^ zq, (p'l / Ae-^'^^T(jX,(y)4A(0)) b) + ^"^"^ (pVZ aWIp) • (1-30) 



The matrix element of the time-ordered product of the two currents is related to the RMC 
amplitude by the time reversal. The form factor gp is contained on the right-hand side 
of Eq. ( |1 . 29| ) in the matrix element of the axial current. If one admits that in the soft 
pion limit only the nucleon Born terms contribute to the divergence of the current-current 
amplitudes, then one has the pion production amplitude which can provide information on 
gp. However, when one of the currents is axial, a contribution to the divergence of the 
current-current amplitudes from the pion pole term in the t-channel survives even in the 
soft pion limit [p^^ -|60[|. A part of this contribution cancels the induced pseudoscalar term in 
the axial current and the remaining part is just the pion pole production amplitude, as one 
can expect intuitively. Then in the soft pion regime, the reaction Eq. (|L29|) is suitable to 
study the weak axial nucleon form factor Fa(/c^) and the electromagnetic form factor (the 
electromagnetic radius) of the charged pion, but not to extract any information on gp. 



In our opinion, the RMC reactions and particularly reactions (|1.23|) and (|1.27|) are at present 



the only available tool to study the form factor gp as a function of the momentum transfer. 
In order to compare our effective form factors with the results of Ref. [0], we define in 



Sect.y the effective Hamiltonian analogously and we consider the velocity independent 



part only. Then in Sect.pl, we present the results for the form factors following from our 
amplitudes [Q up to the order 0{1/M'^). Further we deal with the contribution to gi^s from 
the A excitation amplitudes of our model and we compare our effective weak A^A vertex 
with the one used in Ref. Finally, we discuss the RMC amplitudes stemming from the 
anomalous Lagrangian. 



9 



In Sect.|^, we give the numerical results for the capture rates and we present various photon 
spectra. Our conclusions are presented in Sect.0. Our main result is an enhancement up 
to ~ 7-15 % relative to the calculations without including the A isobar, which can explain 
twice as much of the discrepancy in comparison with the standard approach between the 
gPCAC extracted from the experiment p7| , |28[] , if we take the values of the parameters 

of the model found from the data on pion photoproduction [^^-|47[|. If we are allowed to 



change these parameters independently, we can generate for the value of gp ~ Qp'"^'" the 
photon spectrum, which is in the region of the photon energies > 60 MeV close to the 
spectrum, that should correspond to the experimental one. 

Our triplet capture rate agrees with the one calculated very recently by BHM up to 10%. A 
half of this discrepancy can be attributed to an incorrect integration over the phase volume 
in Ref. pi. 



II. EFFECTIVE HAMILTONIAN FOR RMC 



In presenting the effective Hamiltonian, we follow Rood and Tolhoek |]TP[ . Then the velocity 
independent part is 



H 



(0) 



(1 - (J; • z>) [gi {ai ■ e) + g2{a ■ e) + Qsi^a ■ ex oi] 



+g'A {a'l ■ e) [a ■ kj + g'l [ai ■ e) [a ■ 0) + g'^ i^ai ■ kj {e- 0) + g'^ {e ■ z>) 
^g'^i (^a ■ k X ej + g"i{a ■ u X e) + g'g (0) ■ k^ {a ■ e) + g'l {ai ■ u) {a ■ e) 
+99 {<^i ■ <?) + ^10 ■ ^) ^) + 9io ■ i") '>) 



+g'n ((Ti -k) (a ■ k) {e- v) + g'l^ (ai ■ k) {a ■ z>) {e- 9 



(2.1) 



Here ai {a) are the lepton (nucleon) spin Pauli matrices and z> (k) is the unit vector in the 
direction of the neutrino (photon) momentum vector z7 {k). Not all the form factors are 
independent. Using equation 



iX [k X exj 



1 

V2 



i - \j 



(2.2) 



one gets redefinitions: 

6-2 ^ 5-2 - A {g'j + yg'j) , g[Q 



a'lo + , 5-8 ^ fi'8 + Afl-s , 5-4 5-4 - ■. (2-3) 



where y = {0 ■ k). The last two terms in (|2.1|) are new in comparison with |T^ 
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III. CONTRIBUTION TO H^^p, FROM THE AMPLITUDES OF THE CHIRAL 
LAGRANGIAN OF THE NA-KpujAi SYSTEM 



Here we discuss our amplitudes and contributions to giS. We start by presenting briefly the 
part of the RMC amplitude derived earlier without A's, referring for details to Sect. 3 
of that paper. Then we deal with the amplitudes describing the A excitation processes and 
we compare our effective vertices with those of Ref . . Finally, we discuss the amplitudes 
stemming from the anomalous Lagrangian. 



A. The RMC amplitude without A's 



Besides the muon radiative part M"(/c, q), the amplitude T°'{k, q) 0] consists of three terms 
representing the hadron radiative amplitude 

T\k,q) = ^ {M"(fc,g) + /^(O)e.(A;) [Mf;'^(fc,g) + M;,{7r;k,q) + M^^,{a,;k,q)]} , 



(3.1) 



The amplitude Mj^j)°'{k^q) consists of the nucleon Born terms and of some related contact 
amplitudes. The amplitude M^^^^n; k,q) contains the mesonic amplitude M^-'^ '"(7r; A;, g) 
and all contact terms where the electroweak vertex is connected with the nucleon by the 
pion line. The amplitude M^^(ai; k,q) has graphically a similar structure as the amplitude 
M^^{7t; k,q) with the pion line changed for the Oi meson one. These amplitudes satisfy 
separately continuity equations when contracted with the four momentum transfer of the 
weak vertex. 

Since our model respect vector dominance and PCAC, the sum of the hadron radiative 
amplitudes satisfies exactly the following Ward-Takahashi identities 

g^fM^'^ + M;,(7r) + M^^-,'-^^{a^)] = zUmlAUq) M% + te'^'uip') J'w,M u{p) , (3.2) 
k^ [M^j- + M;,(vr) + M-J-'%a,)] = ie'^\{p') 4,,(gi) u{p) . (3.3) 



Besides, the monopole electroweak form factors with my = nip and tua = appear 
naturally in our amplitudes. Let us note that Eq. (|3.3|) guarantees the gauge invariance of 



the model. The consistency condition |^ for our amplitudes is 

A^(g) k,Ml , = A^(gi) le^^^'M^ . (3.4) 
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Here A^^ ^ is the radiative pion absorption amplitude, is the pseudoscalar vrA^A^ vertex 
and q = k + qi. 



The leading amplitudes are the nucleon Born terms Ml^^"-{k,qY (corresponding to 
M(6), M(c), M{d) in Ref. [|T9l), the amplitudes M;j^(7r,l), M;j^(7r,2) and M^^''*(5) (the 
sum of them corresponds to M(e) in Ref. and the mesonic amplitude M^ '^''" (in cor- 

respondence with M(/) in Ref. [0). As discussed above, besides these amplitudes, other 
contact terms appear. 

The low energy theorems allow one by applying current conservation and PCAC to a 

general amplitude, to determine consistently the amplitude for the RMC in terms of elastic 
weak form factors and pion photoabsorption amplitude, up to terms linear in k and q. As 



shown in Ref. higher order terms cannot be predicted. Since our amplitudes satisfy 
exactly current conservation and PCAC, we can obtain terms of any desired order. We now 
present the expansion of our non-resonant amplitudes up to the order C(l/M^). 



1. Corrections up to the order 0{1/M) 



The non-relativistic reduction of these amplitudes yields the following contributions up to 
the order 0{1/M) to the form factors gi 



9i 

92 
93 

94 

95 
96 

97 

98 
99 

9w 



-K9v 



1 + 



• k 



9vV + 9A>'Vf^ 



N ' 



2M 

-K9a- 9pV + 9v^Vf^v 
-^+9 A + 9mV + 9vV- 9a ^V^^s 



9lV: 



K9a- 9v^V^^v + 9p 



V 



A- + — (1-1/)Ah 

Vflu 



9', 



94 



2M 



9'i 



9i 



2M 



95 



K9v + 9MVf^v + 9v Vf^v + 9m^V, 95 
K9v + i9p + 9l)^Wv + 9mV - ^?y^(l + 2/i 



2M 



V 



K{9v + 9m) + {9p + 9l)^Wv -9vV: 

k 

-K{9v + 9m) - 9v^Vf^s ,98 = 98 
Ki9v + 9m) + i9v + 9M)^Vf^s - 



97 = 97 



2M ' 

'^9AVf^n 



9l 



98 



96 ■ 

k 

2M ' 
u 

2M 



96 



2M 



9'^ 



97 



2M 



9^V 



AMu 



ml + {q 



L\2 



+ K9p- 



m,. 



9io = 9io 



k 
2M 



99 = 9e 



9io 



2M 



9io 



2M 



^For notations see Sect. 3 of Ref. 
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911 = Kg'p^^ , = 9u^ , g'l, = gn^ ■ (3.5) 



Here our notations mostly follow Ref. [p!9f 



771 1 

g = k + u, rj = A± = -(1±A) 

' ' 2M' ^ 2^ ' 



In addition we have 

IJ'V = '^ + f^p-f^n = l + nv, f^s = l + Aip + /is = l + (3.6) 

Besides the obvious momentum dependence of the form factors Qp^^^ given in Eq. ( p..24| ), 
all other nucleon weak vector and axial-vector form factors are assumed to have either the 
monopole momentum dependence, which naturally appears in our model, with my = rup 
and niA = ma^^, or, for the sake of comparison, the dipole one with my = 0.843 GeV and 
ruA = 1.077 GeV B. 



2. Corrections of the order 0{1/M'^ 



Here we have two groups of contributions. The first one arises from the expansion of the 
amplitudes considered above by one order more in 1/M which leads to 

^gi = gM {k - fJ'vi^ ■ k) - (pvl^s - gl^vi^n) (? ■ k 



'2M' 



= -gl\ {v-k)^ 2/x„ [gl + - gf) {s ■ k 

i^gz = -guk + '^gA^i^n {^-k) - '^gyi^n {i^-k 
^g'l = -ghAiJ'S -'^gv^iJ'n, 
^g's = gv^^^ 



"n 1 



-] = (g^-g^)^^n. (3.7) 



M 



The main part of the contribution to the photon spectrum arises from the terms proportional 
to gA and gy in Ag2. These terms appear due to the neutron recoil induced by the time 



13 



component of the weak current. Actually, the terms Ag'l, Ag'l and g'lQ contribute in the 
order 0{1/M^). We have verified that they change the singlet capture rate by ^ 10% and 
the triplet capture rate by ~ 0.8%, which is the reason to keep them. They also arise 
presumably from the neutron recoil. 



The second group of corrections to order (9(l/m^) ^ 0{1/M ) (the hp correction) stems 
from some contact terms present in the hadron radiative part of the amplitude Eq. ( |3.1| ). It 
is discussed in Sect. 4 of Ref. 0. Here we quote the results of the non-relativistic reduction 

As. = -2^^f?^)^„A, = (3.8) 



B. The RMC amplitude with A's 



We derive the RMC amplitudes arising due to the A excitations from chiral Lagrangians 
J^-^. They correspond to the standard nucleon Born terms with the A isobar instead of 



nucleon in the intermediate state. The needed Lagrangian reads 



(3.9) 



Here T is the operator of the transition spin. Another possible term in the pNA vertex is 
suppressed by one order in 1/M and it does not contribute in any sizeable manner ||26[| . 

We take the operator 0^y{B) in a form ||43|-^ 

0M-(5) = 5^, + C(5)7^7,, (3.10) 
C{B) = ]^{l + AB) A + B. (3.11) 



A choice A = —1 simplifies considerably the propagator of the A. 



The coupling constant f-^NA is not well known and the values for f^j^^/An from the interval 
between 0.23 and 0.36 can be found in the literature [Q. From the dispersion theory []61 



/^7VA/47r ^ 0.30 and UnaI^t^ ~ 0-35 from the decay width [0. Also a good fit to the 33 
phase shift was obtained in Refs. ||4^ , |45[| by using f^i^^/An ^ 0.314. The new data on pion 
photoproduction prefer /^tva/^^t ~ 0.371 [0. The ranges of the other relevant parameters 
of the model are |H5|- 



0.8 < Z < 0.7, -1.25 <Y< 1.75, 1.97 < d < 2.65. 



(3.12) 



14 



Our radiative amplitude with the A excitation can be written analogously with the nucleon 
Born term M^j°(l) || as 



M,^;'^ = -u{p') 



u{p) , 



(3.13) 



Here the weak A^A vertex reads 



(3.14) 



with the vector part defined as 



and with the axial-vector part of the form 

f-n-firNA 



'JA,fiaiq) 



miA;\{q)-q,qxAl{q) 0^x{Z) 



(3.15) 



(3.16) 



Further the electromagnetic 7A^A vertex is 

Jem,u'y{k) = -Jv,,yy{k, /c^ = 0) = -i (j^^ {kflS^X - kx^i,^) 0^(j{Y) '^r^lx ■ (3.17) 

Finally, S'P {p) is the A isobar propagator. With the choice Oa/3 = dap, our amplitudes 
Eqs. ( 3.14| )-( |3.16| ) coincide in the form with the ones obtained in Ref. from the study 
of the weak A^-A vertex in the reaction v d ^ A++ n [|. 



For the divergence of the resonant amplitude M^j"" from Eq. ( p.l3| ) we have 



(3.18) 



Here the associated resonant radiative pion absorption amplitude M^'^ is 



M. 



A, a 



-u{p') 



uip), (3.19) 



^For a recent study of this reaction see |63|. 
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and the vrA^A vertex reads 



(3.20) 



We now present the contributions from the amphtude Eq. (|3.13|) to the form factors Qi. They 
are 

A^i = \kC- - C^)C {1 + (1 - R) [CiY) + C{Z) + 2(2 + R)C{Y)C{Z)]} , 
^92 = \{C+ + C_)C(1 - R) [-(1 + 2R) + 2(1 - 2R)C{Y) + 2(1 - R)C{Z) 



+A{2- R)C{Y)C{Z)] + {C+ + C^)C 



9^ 



6Mg 



N 



{-[{1- R)il + 2R)k 



+yu] + 2(1 - R) {[(1 - 2R)k + yu] C{Y) + [(1 - R)k + yu] C{Z) 
+2 [(2 - R)k + (2 + R)yu] C(r)C(Z)}) , 

^93 = ^A(C+ + C_)C{l + (l-i?)[C(r) + C(Z) + 2(2 + i?)C(r)C(Z)]} , 
A^?; = -Ag', = (C+ + C_)C, 



^3M<?^ 



(C+ - C_)C {1 + (1 - i?) [C(F) + C(Z) + 2(2 + R)CiY)CiZ)]} , 



Ag[, = (C+ + C_)C 



6M< 



{1 - 2(1 - R) [C{Y) + C(Z) + 2(2 + R)C{Y)C{Z)]} , (3.21) 



where 



C 



4 fnfTTNA 

3 



(Ma - M) + 
(Ma - M) + 



2M 



M/M^ 



Ma + M 

2M 
Ma + M 



m, 



m,, — V + 
^ 2M 



2v — — 



2i^ 



+ yA;) 



z/ - + 



2M 



(2z/ - m^) 



(3.22) 



According to the concept developed in Refs. [43-45|, we take the mass of the A isobar real. 



C. The RMC amplitude from an anomalous Lagrangian of the npujai system 



We have considered so far the amplitudes where a natural parity does not change in any 
vertex. The natural parity of a particle is defined as P (— 1)*^, where P is the intrinsic parity 
and J is the spin of the particle. Some amplitudes of this kind relevant for the process under 
study are presented in Fig. |I]. 
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FIGURES 



P' 



P 




P' 



P 




P' 



P 




FIG. 1. (a),(b) - the radiative hadron amplitudes obtained from the anomalous Lagrangian of 
the npujai system, Eg. ( fj.23 ); in (a), B=7r or ai; (c) ~ the associated radiative pion absorption 
amplitude. These amplitudes satisfy PC AC. 



The starting point is an anomalous Lagrangian of the npuai system constructed 
within the approach of hidden local symmetries |p3JBB] . The electromagnetic interaction in 
such a system was first considered in Ref. and the relevant constants Cj were extracted 
from the data as well. The weak interaction was incorporated explicitly in Ref. and the 
refit of the constants to the modern data was made in Ref. |^ . 



The Lagrangian reads 

( 1 



(3.23) 



where besides the meson fields, the external vector isoscalar B^ and isovector and axial 



vector isovector fields are also included. The constants q are |^| 

C7 = 8.64 X 10"^ , Cg = -L02 x 10"^ , cg = 9.23 x 10"^ , cio = 1.29 x 10"^ . (3.24) 



The axial RMC amplitudes arising from the anomalous Lagrangian Eq. (|3.23| ), are 
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2 



2M^ 



(3.25) 



and they correspond to the processes presented in Figs.|I|(a) and 0(b). Together with the 
radiative pion absorption amphtude of Fig.|l|(c) 



^ana 

TT, V 



9p 



(3.26) 



the amphtudes Eq. ( ^1251) satisfy PCAC: 



5m 



M-''^(l) + M-''^(2) = iUmi/\Uq)M:y 



(3.27) 



The contribution from the anomalous amphtudes Eq. ( |3.2 
amphtude Eq. (|3.1|) for the reaction Eq. (|1.23| ) is given as 



to the total hadron radiative 



(C7 + Cg) 



(1) + m;:''^ 



(2)' 



(3.28) 



In Ref. [Q, a contribution arising from the Wess-Zumino-Witten anomalous Lagrangian 
was estimated. Graphically, it corresponds to our Fig.|l|(b) with the pion instead of the rho 
meson and with the vector interaction instead of the axial one. The associated amplitude 
depends on the momentum transfer q^ and therefore, it does not possess the enhancement 
factor ^ 3 for large photon momentums. For illustration, we present the contribution to 
one of the form factors 

Ag', = -n ^lJ'\-^y% ',, (3.29) 
87rV^ 2m^ 



all other contributions have a similar structure. It is also seen that these contributions are 
of the order 0{l/M^), because M'^fl ~ M^. 

In our case, it is the amplitude related to the graph Fig.|l|(a), which is g^-dependent. We 
present from the calculated contributions to Qi arising from the amplitudes Eq. ( p.25|) only 
the one for the form factor (72, the other ones are suppressed by one order in 1/M 

gl f2MY, , k ^ gl (2mV ^ , k-s 



18 



For comparison, we keep also the (?p -dependent contribution. Using the KSFR relation 

2 



2f^gp = m?, we can rewrite this contribution to the form 



1+ Ky k 



9 A 'M12/2 



Taking into account that 12/^ ^ M^/10 we can see that the ^fp-dependent contribution 
is larger than the (yfp-dependent one for large values of /c by a factor ^ 20 . However, it is 
not enough to influence the photon spectrum, because of an additional factor (cy + cg) in 
the amplitude Eq. ( p.28| ) (see below). It is seen from Eq. ( p.30|) that the first term on the 



right-hand side is suppressed in comparison with the second one arising from the contact 
amplitude M^"'"(2). As we shall see later, the second term contributes to the triplet capture 
rate by an amount of ~ —0.2%. 

Let us note that the sum of the vector RMC amplitudes arising from the anomalous La- 
grangian Eq. (|3.23|) is zero with a good accuracy. 



One can obtain more general result for A5f2, Eq. ( |3.30D , by a change 

1/3 9ty\l9w\^ i^v 9p2/9pi, 



which coresponds to the model used in Refs. [44,45] for describing the pion photoproduction 



amplitude in the t channel. In this model, the 7rp7 and ttcvj amplitudes are effectively the 
same as those obtained from our anomalous Lagrangian Eq. ( |3.23| ) and the pNN and uNN 
vertices contain four constants gpi, gp2, g^ji and g^^21 which are the free parameters, obtained 
together with other free parameters of the model from a fit to the data. 

Compared with the sets of the form factors given in Eqs. ( ^.7] ) and (|3.8|), the giS, Eq. ( ^.3U| ), 
are even larger. However, due to the values of Cj, Eq. (|3.24|) , the factor (C7 + C9) ~ 1.8 x 10~^ 



makes the contribution from the amplitude T"„, Eq. ( |3.28| ), small. 



IV. RESULTS 



Using the Hamiltonian -f^e//' ■'^l- (|2-1|) , and the sets of the form factors (^j, Eqs. (|3.5| ) 



(^.7D,(P78D, and ( |3.21D , we have calculated the capture rates and the photon energy spectra 
for the RMC in a muon-hydrogen system described by a spin density matrix = s, t) 



dk 



— (a'^Gp cos 9cm/mn) mMnk [ dy 



47r 



W + k{y-l) 



Tr{(l - a- z>) 



H. 



(0)- 

eff 



(4.1) 



19 



Here a is the fine structure constant, Gp\^ the Fermi constant, cos 6c is the Cabibbo angle, 
m is the reduced mass of the ^ — p system, the neutrino energy is determined by the energy 
conservation 



W 

^0 = W + k{y- 1) ^^'""^ " 



X (m^ + k{y - 1))] (A;„„^ - k) , (4.2) 



where the maximum photon energy is given as 



2W 



kma. = " ' W = Mp + TTl^, (4.3) 



and Mp(„) is the proton (neutron) mass. The singlet and triplet spin density matrices are 
3|,|68| 



Ps = \{l - S-Bi), = 1(1 + lc?-aO. (4.4) 



We have also calculated the capture rates and spectra for the ortho- and paramolecular p/xp 
states and for the mixture of muonic states relevant to the TRIUMF experiment p7| , p8| . 
The ortho (Ao)- and paramolecular (Ap) capture rates are given in terms of A^ and At as 



Ao = 0.756 A, + 0.253 At, Ap = 0.286 A, + 0.857 At , (4.5) 



and the capture rate Ay, relevant to the TRIUMF experiment ||27| , |28 



IS 



At = 0.061 A, + 0.854 Ao + 0.085 Ap. (4.6) 



Now we present numerical results for the capture rates. 



A. Capture rates 



Here we present the results for the capture rates calculated in various models. If not stated 
otherwise, we use the monopole form factors and we put x = 1 in Eq.( 1.24 ). 



A) We first discuss the results obtained in the model with the A isobar kept on-shell. We 
give the singlet and triplet capture rates in a more detail, in order to see explicitly various 
contributions. 
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A, X 10^ = 0.40(0) + 1.65(-1) + 1.29(-2) + 0.11{hp) - 0.07(a/) + 0.05(A) 

= 3.43 (3.51) (4.7) 

At X 10^ = 43.7(0) + 53.1(-1) + 3.7(-2) + 0.9{hp) - 0.2 (a/) + 2.2(A) 

= 103.0 (103.4) (4.8) 



Here on the right-hand side of Eqs. ( [4.7|) and (|4.^ ), the number n (n=0,-l,-2) in the brackets 
means the order of the contribution in (9(1/M") and hp (al) and A mean the contributions 
from the hard pion form factors Eq. ( p.8|) (from the form factors Eq. ( p.30| )) and from the 
form factors Eq. ( |3.21| ) due to the A isobar excitation processes, respectively. These were 
calculated using the parameters 



2 



4:71 



0.371, Gi = 2.525, Y = Z = -0.5 . (4.9) 



The choice of the parameters Y and Z is such that only the terms proportional to 5^1, in 
Eq. ( |1.28|) contribute (the A is on-shell). The contribution of the A excitation to is 



«i 2%, which is in agreement with BHM. The numbers in the brackets on the right-hand 
side of Eqs. (^4.7|) and (|4.8|) are obtained using Eq. ( |1.24| ) for gp^^^ without the correction 
^p*-^-* included. 

For the other capture rates we have 

Ao = 28.7 X lO^^s"^ , Ap = 89.3 x lO^^s"^ , A^ = 32.3 x lO-^s"^ . (4.10) 



For the dipole form factors, analogously with Eqs. (|4.7| ) and (|4.8| ) we obtain 

A, = 3.28 (3.50) X 10" , At = 101.5 (102.5) x 10"^s"^ . (4.11) 



Let us compare our results with available calculations. The singlet and triplet capture 
rates were calculated earlier by Opat [1^. He obtained A^ = 4.96 x 10~^s~^ and At = 
90.0 X 10-V^ 



Very recent calculations yield A, = (2.90-3.10) x lO^^s"^ and A* = (112-114) x lO^^s^^ 
Having in mind that A^ results as the difference of two large and almost equal numbers, the 
agreement between our value Eq. ( [4.7|) of A^ and the BHM value can be considered as 
satisfactory. However, the difference of ~ 10% between the triplet capture rates is too large. 
A half of this discrepancy can be understood by checking the integration over the phase 
volume. We use for the neutrino momentum Eq. ( [4. 21) with k^ax from Eq. ( [4. 3]) , while BHM 
employ for kmax equation (4.37), which is in our notations 

. ^ , ^)-. . - ^ . ^) . (4,12) 
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where Mn = {Mp + M„)/2 is the nucleon mass. From Eg. ( |4.3|) , one obtains k^ax = 
99.15 MeV, while from Eq. ( ^l2D one has k^ax = 100.3 MeV and instead of At = 103.0 x 
10~^s~^ one obtains = 108.0 x 10"^s"^, which is larger by ^ 5%. Using in Eq. ( ^^2] ) for 
uq the expansion Eq. ( [4.12|) for kmax, one obtains 



I'o = m^-k 



m, 



+ 



k 



X 



2Mn Mn 
mi + 2k{m^ - k){y - 1) 



(1 - i/)(m^ -k) + 



(4.13) 



In Ref. Eq. (4.39) is used for z/q. It retains terms up to the order 0(l/MAr)0, which 

yields in our case Aj = 107.2 x 10~'^s~^. But the source of the remaining difference of ~ 5% 
between the results for the A^ is not clear. 



B) The results of calculations without the A excitation effect are 

A, = 3.38 X lO^^s"^ , kt = 100.8 x lO^^s"^ , At = 31.6 x lO^^s"^ . (4.14) 



C) We now present the capture rates for the same case as in B), but for the value of the 
parameter x = 1.5. The rates are 



A, = 8.35 X 10"^ , At = 116.6 x 10"^ s'' , At = 39.8 x 10~^ . 



3 „-l 



-3 „-l 



(4.15) 



The strong dependence of the capture rates on gp was already known to Opat |18 



D) The capture rates, calculated for the same parameters as in A), but Y = 1.75 and 
Z = —0.8, are 



A, = 3.28 X 10"^ s'^ , At = 106.0 x 10~^ s'' , At = 33.0 x 10"^ s 



3 „-l 



-3 „-l 



(4.16) 



In this case, the effect of the A excitation is for the values of Y and Z, allowed by the 
inequalities of Eq. ( p.l2|) , maximal and it is ~ 5% for A^ and ^ 4% for At. 



E) The capture rates, calculated as in D, but Z = —1.95, 



-3 „-i 



= 5.93 X 10"^ , At = 116.7 x 10"'' s'' , At = 38.1 x 10"^ s 



(4.17) 



These capture rates are close to those calculated in the case C) for x = 1.5. In order to 
achieve the same enhancement in the rates, one should put A isobar more off-shell than it 



There is factor 2 missing in the denominator of the third term on the right-hand side of this 
equation. 
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is needed in pion photoproduction processes. On the other hand, the values of the off-shell 
parameters X, Y, Z depend strongly on whether the pion production amplitude is unitarized 
or not and on the method of unitarization. Let us note that our choice of the parameters of 
the model is not optimal. In order to extract an optimal set of these parameters from the 
data, one should use a minimization procedure. 

Similar enhancement in the capture rates as in the case E) can be achieved by taking 
f^N/^/4:7r ^ 20 or Gi ~ 20, which is an amplification of ~ 7 in the ttA^A coupling and of 
^ 8 in the constant Gi, which is much more than the enhancement factor of ~ 2.5 for the 
parameter Z. 



B. Photon spectra 



Photon spectra, corresponding to the capture rates calculated in the model A of the previous 
section, are presented in Fig.|^. The A isobar is included, but the choice of the parameters 
Y = Z = —0.5 is such that the isobar is on-shell. As it is seen from Fig. ^ our spectra are 
in a close correspondence with those of Fig. 3 of Ref . p6[ . However, our spectrum for the 



triplet state of the fi—p system (long-dashed curve) differs from the analogous spectrum of 
Fig. 5 as it should be, because the triplet capture rates differ significantly. 
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FIG. 2. Photon spectra calculated in the model A of the previous section, the A isobar 
excitation effect is included; the dotted, long-dashed, dot-dashed and dashed curves correspond, 
respectively, to the singlet, triplet, ortho- and para -piip molecule spin combinations, the solid curve 
correspond to the mixture of muonic states relevant to the TRIUMF experiment [27,2^]. 
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FIG. 3. Relative difference in % for the spectra calculated with (model A) and without (model 
B) the A excitation included; the dotted, long-dashed, dot-dashed and dashed curves correspond, 
respectively, to the singlet, triplet, ortho- and paia pfip molecule spin combinations, the solid curve 



correspond to the mixture of muonic states relevant to the TRIUMF experiment [27,2 
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FIG. 4. Relative difference in % for the spectra calculated with (model D) and without (model 
B) the A excitation included; the dotted, long-dashed, dot-dashed and dashed curves correspond, 
respectively, to the singlet, triplet, ortho- and parap/ip molecule spin combinations, the solid curve 



correspond to the mixture of muonic states relevant to the TRIUMF experiment [27,2 



The percentage change in the spectra when the A excitation effect is taken into account is 
presented in Fig. ^ The A excitation effect was calculated according to the model A of the 
previous section. The case without this effect corresponds to the model B. This change in 
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the spectra due to the A was first calculated by Beder and Fearing p6[. Our Fig.|^ is in a 
good agreement with Fig. 4 of Ref. |26], but it differs from the analogous Fig. 6 of Ref. [29 . 



Similar calculations are presented in Fig. ^ and Fig. |^ using instead of the model A the 
spectra of the model D and E, respectively. As it is seen, by putting the A isobar off-shell, 
both the singlet and triplet spectra are changed sensibly. 
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FIG. 5. Relative difference in % for the spectra calculated with (model E) and without (model 
B) the A excitation included; the dotted, long-dashed, dot-dashed and dashed curves correspond, 
respectively, to the singlet, triplet, ortho- and para pfip molecule spin combinations, the solid curve 
correspond to the mixture of muonic states relevant to the TRIUMF experiment 
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In Fig. ^ we show how the photon spectra, relevant to the mixture of the muonic states 
for the TRIUMF experiment, depend on the parameters of our model. The dotted curve 
corresponds to the model B (no A, x=l), the dashed curve is the photon spectrum for the 
case A (A on-sheh, Y = Z = —0.5), the dot-dashed curve is calculated using the model 
D: A is off-shell, the parameters Y = 1.75, Z = —0.8 are at the boundary of the region 
Eq. ( p.l2|) allowed by the pion photoproduction data ||4^-^. In this case, about two times 
more discrepancy is explained in comparison with the dashed curve. The dependence of 
the photon spectrum on the change in gp is illustrated by the long-dashed curve, calculated 
within the model C (no A, x=1.5). We have found that this curve, normalized to 14.5 counts 
at /c = 60 MeV, follows closely the solid curve of Fig. 4 PB| . All other curves in this and in 
the next figure are multiplied by this normalization factor. The long-dashed curve closely 
follows the solid curve from A; ^ 60 MeV, which corresponds to the model E: A is off-shell, 
the parameter Y = 1.75 is the same as in the case D, the parameter Z = —1.95. 



Evidently, this picture is in agreement with the one obtained by comparing the Beder- 
Fearing model |2^ with the experiment |2^j2^: the experimental photon spectrum can be 
satisfactorily described with gp of the form Eq. ( |1.24| ) only for a; ~ 1.5 and only a small part 
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of the discrepancy could be explained by calculating the spectra by including the on-shell 
A isobar. As it is seen from the solid curve, putting the A isobar off-shell, one can describe 
the experimental data quite well. 




As we have just found, our model can explain about two times more of the discrepancy 
between the gp'~^^'~^, Eq. ( |1.14| ), and of gp extracted from the experiment p7| , |28| , if one 



restricts oneself to the set of the parameters from Eq. (|3.12|) . They were extracted in Refs 



^^-^T] from the data on pion photoproduction by using unitarized multipoles arising from 



the pion photoproduction amplitude. The problem with the unitarity appears in the pion 
photoproduction, because of the pion-nucleon interaction in the final state. It is true that the 
time reversed of the pion production amplitude is connected with the hadron radiative part 
of our RMC amplitude by the continuity equations Eqs. ( p.2|) , (|3.18|) , and ( p.27|) . However, 



we do not need to unitarize our amplitude. Therefore it is not clear, how the inequalities 
Eq. ( |3.12|) are restrictive for the problem considered here. It is seen that the model can 
reproduce reasonably the experimental photon spectrum [0,^ in the region A; > 60 MeV 
for the PCAC value of gp, Eq. ( |1.14|) , if the values of the parameters Y and Z are taken to 
be y = 1.75 and Z = -1.95. 

In Fig.[^, we show the dependence of the calculations on uncertainties in our knowledge of 
gp and of the admixture ^ of the 5 = 3/2 orthomolecular pfip state. As discussed recently 
||33|| , admixture of the 5 = 3/2 orthomolecular pfip state changes the molecular capture rate 
to @ 

a:, = eAo(l/2) + (l-OAo(3/2), (4.18) 
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where Ao(l/2) is of Eq. (U) and Ao(3/2) = 1.009 At. It was found in H that the data 
on OMC in hydrogen requires gp < 1.2 ^fp*^^*^ or ^ > 0.95. This restriction on gp is in 
agreement with our Eq. ( |1.20| ) for the OMC in ^He. The dependence on the uncertainty in 
gp and on ^ is illustrated by the long-dashed, dot-dashed and dashed curves. The numbers 
in the brackets are for the total capture rate, for the partial capture rate in the interval 
(60 — kmax) MeV and for the capture rate in counts for this interval, respectively. Otherwise, 
the unit for the capture rates is 10~^s~^. For the sake of illustration, we assigned a 10% error 
p8| to a set of 'data' represented by 14 points of the spectrum C. The number of counts 
276, related to this curve is to be compared with the number of counts 286, which can be 
read off the histogram presented in Fig. 4 ||2^ . As it is seen, all curves lay already inside the 
2cr bound. 
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FIG. 7. Dependence of the photon spectra, corresponding to the mixture of muonic states 
relevant to the TRIUMF experiment [ p7| , p^ ], on various parameters of the calculations; the dotted 
curve is for the reference and corresponds to the model A (32.3/7.1/196); the long-dashed curve is 
calculated for the model D, but ^ = 0.95 (36.3/8.1/225); the dot-dashed curve is obtained within 
the model D, but x = 1.2 (35.9/8.6/239); the dashed curve- the model D, but x = 1.2 and ^ = 0.95 
(39.4/9.3/259); the solid curve corresponds to the model C (39.7/9.9/276). For the details see the 
text. 



V. SUMMARY 



In this paper, we have presented the capture rates and the photon energy spectra for the 
RMC in hydrogen, calculated using the effective Hamiltonian -f^e//' -^<1- ( P-IP ? where the 
form factors gi are obtained from the amplitudes derived from the chiral Lagrangian of 
the NAnpujai system. The non-resonant part of the Lagrangian contains the normal and 
anomalous Lagrangians of the Nnpujai system interacting with the external electromagnetic 
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and weak fields by the associated one-body currents Pfl ^ , p2| , |53| . In the expansion of the 
amplitudes in 1/M, we keep all terms up to the order 0(1/M^). 

For the resonant part of our Lagrangian, we have extended the standardly used model 
P J26|J29|J4^ by adopting results of the model developed by Olsson and Osypowski |^ and 



by Davidson, Mukhopadhyay and Wittman |4^,^, allowing the A isobar to be off-shell. 

The calculated triplet capture rate differs by ~ 10% from the one derived quite recently 
within the HBChPT approach in Ref. pUf. About a half of this discrepancy can be under- 



stood by the use of an approximate equation for the neutrino momentum in Ref. [^]. The 
origin of the rest of the discrepancy is not clear. 

In the model, restricting the A isobar on-shell, our spectra are close to those obtained earlier 
by Beder and Fearing Our full model, that includes off-shell A isobar, can explain two 
times more of the discrepancy between the PC AC value of gp, Eq. ( |1.14|) , and of gp extracted 



from the experiment [^,^], if one restricts oneself to the values of the parameters of the 
model extracted from the data on the pion photoproduction off the nucleon. Let us note that 
taking into account existing uncertainty in g'p'^^'^ extracted from the OMC in the hydrogen 
and ^He and in the parameter ^, regulating the admixture of the S = 3/2 orthomolecular 
pfip state, one finds that this model is, actually, in reasonable agreement with the data 
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Moreover, if one is allowed to vary the parameters of the model independently, the exper- 
imental photon spectrum can be described for the induced pseudoscalar form factor gp of 
Eq. ( |1.24D without any scaling (x = 1). It would be difficult to find any physics behind the 
scaling of gp, which would mean violation of PCAC. On the other hand, the part of our 
full model, containing the electroweak interaction of the off-shell A isobar, is widely used 
to describe successfully the pion photoproduction data, which provides a strong basis for 
confidence in our results. For the reaction of the RMC in the hydrogen, Eq. it is the 



only model known so far, providing enough enhancement in the high energy region of the 
photon spectrum. 

In conclusion we note that the reactions of the RMC in the hydrogen and ^He are at present 
the only available effective tool for the study of the form factor gp as a function of the 
momentum transfer. Therefore, more efforts, both theoretical and experimental, are highly 
desirable. 
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